Alfven modes driven nonlinearly by metric perturbations in 
Anisotropic Magnetized Cosmologies 
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We consider anisotropic magnetized cosmologies filled with conductive plasma fluid and study 
the implications of metric perturbations that propagate parallel to the ambient magnetic fleld. It 
is known that in the first order (linear) approximation with respect to the amplitude of the pertur- 
bations no electric field and density perturbations arise. However when we consider the non-linear 
coupling of the metric perturbations with their temporal derivatives, certain classes of solutions 
can induce steeply increasing in time, electric field perturbations. This is verified both numerically 
and analytically. The source of these perturbations can be either high-frequency quantum vacuum 
fiuctuations, driven by the cosmological pump field, in the early stages of the evolution of the Uni- 
verse, or astrophysical processes, or a non-linear isotropization process, of an initially anisotropic 
cosmological spacetime. 
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I. INTRODUCTION 

Magnetic fields are known to have a widespread pres- 
ence in our Universe, being a common property of the in- 
tergalactic medium in galaxy clusters [H, while, reports 
on Faraday rotation imply significant magnetic fields in 
condensations at high redshifts [2] . Studies of large-scale 
magnetic fields and their potential implications for the 
formation and the evolution of the observed structures, 
have been the subject of continuous investigation (see e.g. 

for a representative though incomplete list). Mag- 
netic fields observed in galaxies and galaxy clusters are 
in energy equipartition with the gas and the cosmic rays 
[r3 |. The origin of these fields, whether of astrophysical 
or cosmological origin, remains an unresolved issue. 

If magnetism has a cosmological origin, as observa- 
tions of /iG fields in galaxy clusters and high-redshift 
protogalaxies seem to suggest, it could have affected the 
evolution of the Universe [IJI- There are several scenar- 
ios for the generation of primordial magnetic fields (see 
e.g. [IBl)- Most of the early treatments were Newto- 
nian, while relativistic treatments appeared recently in 
the literature. A common factor in all these approaches 
is the MHD approximation, namely the assumption that 
the magnetic field lines are effectively frozen in an in- 
finitely conducting cosmic medium (i.e., of zero resistiv- 
ity). With few exceptions (e.g. [ll] , [l^l , [Hi ) the role 
of non-zero resistivity and kinetic viscosity have been ig- 
nored, these features however being essential for a com- 
prehensive picture of the non-linear magnetic field evolu- 
tion. The electric fields associated with the resistivity can 
be the source for particle acceleration, while the induced 
non-linear currents may react back upon the magnetic 
field [H]. 

Many recent studies use a Newtonian or a FRW cosmo- 
logical model to represent the Universe and super-impose 
a large-scale ordered magnetic field. The magnetic field 



is assumed to be too weak to destroy the FRW isotropy 
and the anisotropy, induced by it, is treated as a per- 
turbation 0, [nj, [l^. Current observations provide a 
strong motivation for the adoption of a FRW model, but 
the uncertainities of the cosmological Standard Model are 
several. Therefore the approximation of neglecting the 
background anisotropy and magnetic fields, may lead to 
effects and phenomena that are absent in the above treat- 
ments. Within this context the formation of small-scale 
structures and the excitation of resistive instabilities in 
Bianchi-type models has been explored long ago [lB| , but 
the issue of excitation of MHD modes in anisotropic cos- 
mological models and their subsequent temporal evolu- 
tion is far from being exhausted [20| . 

Usually one assumes that h? is small. However one 
may consider strong GWs either in the Early Universe 
when dynamical isotropization of an anisotopic cosmolog- 
ical model takes place, or due to amplification of quantum 
vacuum fluctuations from the the variable gravitational 
field (cosmological pump field) of the Universe (metric 
tensor) . Also hh can be large due to high frequency grav- 
itational wave perturbations in the Early Universe. 

On the other hand examples of electrovacuum cylindi- 
cally symmetric spacetimes, interacting with GWs, in the 
hiW theory are well known (^22*1, [23^), or the extraction of 
GWs by black hole collisions [24] whereas the covariant 
decomposition into scalar, vector and tensor perturba- 
tions is treated in (25j . Gravitational waves can carry a 
large amount of energy near the sources where they are 
generated (see e.g. [2g]). Though they do not interact 
much with matter under normal conditions, in the lin- 
ear level, it has been shown that they can excite various 
kinds of plasma waves, more efficiently with increasing 
background magnetic field (see e.g. [2l[, [23|, HI], [1^, 
Isol. I3ll[. [S^l . fssjl. Non-hnear effects have been studied 
i7[23, [3| 

However the problem of the behaviour of conductive 
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plasmas, in anisotropic magnetised cosmologies, which 
are driven by external perturbations, in the full non- 
linear level is more or less open. Motivated by this fact we 
consider this problem in the background of Thome's class 
of anisotropic magnetized cosmologies. We write the per- 
turbed Einstein's equations in a manner that takes into 
account the non-linear coupling of the metric perturba- 
tions with their temporal derivatives. Then we study 
mainly the electric field perturbations induced by classes 
of solutions, to the evolution equations of the metric per- 
turbations (which however are written in the linear ap- 
proximation). The organization of this paper is as fol- 
lows: In Section II we present the derivation of the closed 
set of equations used, aided by the two appendices. In 
Section III we present the various classes of analytical so- 
lutions that enter in our discussion, classified into small 
and large-t solutions. In Section IV we present the nu- 
merical results of our paper and finally in Section V we 
present a brief discussion of the results. 



II. BASIC EQUATIONS 

We consider classes of anisotropic magnetized cosmolo- 
gies with mater content in the form of a perfect fluid and 
with the magnetic field in the z-direction. The metric is 
taken to be 
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where A = A{t), W = W{t)^ for the metric compo- 
nents of the background spacetime (see Appendix A for 
the class of Thome's anisotropic magnetized cosmologies 
that is used), while the metric perturbations are assumed 
of the form /i+ — h^{t, z) and hx ^ hxit, z). 

The energy-momentum tensor for the perfect fluid is 
taken as r^{'^ = T^u' + (^T)i,y) where of course in this 
expansion we use the background value 
T^"J = diag{p{t),A^p{t),A^p{t),W^{t)p{t)) and the per- 
turbation is computed as 



Here we have Sp = 5p{t, z), Sp — Sp{t, z), for the diagonal 
terms, 6To^ ^ -{p + p)A'^{Su''), SToy = -{p + p)A'^{Suy) 
and STqz — —{p + p)W'^{6u^). Also we have used the 
four-velocity perturbations around the background (co- 
moving) value u^Q^ = (1,0,0,0) and perturbation of the 
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structed from the field tensor 
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where E^ = F3^'u^, B'' = {l/2)e^^'"^UbFcd are the elec- 
tric and magnetic field respectively, and the background 
magnetic field is assumed (as in Appendix A) to be 
Ba{t) — BQ{t)z. In general we assume the perturbations 
SE^ = SEi{t,z) and SB^ = SB^t,z), {j = x,y,z) and 
the energy-momentum tensor is 



T, 



(em) 

^p2 .^^ F'^^F^'^g^pg^,) 



(4) 



For the metric of Eq. ([T]) the components of the Ein- 
stein tensor are computed. Then these are expanded up 
to the second order and we set (G = 1 = c) 



(5) 



Here on the r.h.s we have the total energy-momentum 
tensor, both form the fiuid and the EM field, expanded 
with respect to the fiuid and EM-field perturbations. The 
(^G^i/— term contains in general first and second order 
corrections, with resect ot the GW amplitude, (though 
here the first order terms vanish as it is well known for 
the Alfven modes [2l[) and we set, as indicated below, 
SGfj^jy — SnST^i,. Also equating the zeroth order terms we 
obtain the usual Einstein's equation for the background 
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8np + AHBo)^ 
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= %^p-A\B^f (6) 



It is easy to show that the class of cosmological models 
of Appenix A, satisfies identically Eqs. ([6]). Now since 
we have ^(G^S -I- (SG^^) = -f 0{h^) we will also have 

from Eq. ^ ViT^"J + 6T^^) = + 0{h^). So these 
conservation equations need not be considered separately, 
because they are embodied in Eqs. ^ and in 




SG., = 8nST, 



(7) 



The energy-momntum tensor for the EM-field is con- 



The content of Eqs. ([7]) determines the fiuid and EM- 
field perturbations in terms of the propagating GW and 
will be given below. 

The only remaining set of equations that need to be 
considered is Maxwell's equations V^i^'^" — AttJ^ and 
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P[a(3;-y] = 0- Here the current density is given by 
J'^ = pqw^ + aF'^^Ujy. The local charge density pq — 
ZgUi — eue ~ 0, i.e., can be taken as zero due to the 
mobility of the lighter electron with respect to the other 
ion species, or due to exact cancelation (i.e., in the case 
for example of an electron-positron plasma). Also we 
have a = l/(47r?7) for the finite conductivity of the plasma 
fluid, with 7] the resistivity. 

Now from the {aj3; 7) — (Oxy) and (a;?;z)-components 
of the sourceless Maxwell's equations we exactly obtain 
Eq. dSJ) of Appendix A, and 



SB^it,z) = ^ihl + hl), 



(8) 



while from the (Oxz) and (Oyz)— components we obtain 
respectively 

{A''W^6B^),t-A\SEy),, = (9) 

From the sourcefuU Maxwell's equations we obtain that 
SE"" = and 
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along with a similar equation, that results from the sub- 
stitutions 5E'' SEy, SBy -SB"" and SuV -Su'=. 
These are the propagation equations for the EM-field per- 
turbations. 

We now consider Eqs. ([T])- From the (Ox), (Oy) and 
(Oz) components we obtain the fluid's four-velocity per- 
turbations 
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From proper combinations of the (xx), {yy) and {zz) 
components of Eqs. ([7]) we obtain the pressure perturba- 
tion 
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[h+h+ + hy.hy] - ^^[/^ + hi] = 



(15) 



and the propagation equation for the cross polarization 
of the gravitational wave 
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(12) Finally from the {xy) component we obtain the propaga- 
tion equation for the times polarization of the GW 
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□/ix = [8np + A'^{Boy]hx- Eqs. constitute a 

closed system of equations that we study below. 



III. ANALYTICAL SOLUTIONS 

Substituting the solution of Appendix A into Eq. (fT51) 
we obtain (setting also h^{t,z) — e^^^^h^{t)) 



2j \ h 



27 



(17) 



A similar equation is satisfied by the other polarization. 
Also the constraint equation, Eq. (|15p is satisfied if both 
polarizations obey 
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We can take as final time of our analysis any time be- 
fore the recombination epoch tfin — where pre- 
sumably the cosmological plasma ceases to exist. Also 
the initial time can be any time after the typical time 
that inflation starts {to ~ 10~^'*sec), or the end of the 
reheating period {to ~ 10^^^ sec). Between these we rec- 
ognize two eras: The first is the smaZ/-time era where in 
Eq. (|16p the second term is small compared with the next 
two. The second is the large-time era where the last term 
can be omitted. Equating these we have an estimate of 
this transition time as 
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(19) 



We observe that this time scale coincides with the time 
that the mode, with comoving wavevector kg = 1/Xc, 
enters the Horizon. Indeed the physical wavelength 
Xphys ~ WXc must be smaller than the Horizon scale, 
Xphys < Ih — c/Hw = cW/W. This gives the estimate 



[^^/fcg]^/^*^^ ^\ which is of the same order with 



t 

Eq. (HID. 

We identify the modes by their wavenumber kg and 
their initial-time frequency /o so that 27r/o — kg, (we 
take throughout c = 1). At the end of the time interval 
tfin let the frequency be //m = fm{.KHz). Due to the 
cosmological redshift we have /o = //mT^(t/m)/W^(io)- 
If we denote u := (t/to) then the transition time of 
Eq. is given by 
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(20) 



and is plotted in Fig. It is evident that this can be 
obtained for other boundary values for to, tfin. 
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FIG. 1; The transition time t* 
the parameter 7 



of Eq. (|19p . as a function of 



(18) 



A. Small-t Solutions 

There are many types of solutions, for example power- 
law solutions, namely when the second term on the l.h.s 
of Eq. ini) is omitted. The first is 



(21) 



and it is evident that by substituting into Eq. p7|) the 
error is of order ©(l/t^'). We also find that Eq. ^ 
is satisfied with remaining error terms of order 0(l/i^'), 
provided that Ci = 0. 

The second solution we found is 
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which, when substituted into Eq. (|17p gives error terms 
of order 0(1/^^)! Also the constraint equation Eq. 
is satisfied to the order 0(l/i'+^). Here the mode is 
outside the Horizon and the amplitude of the GW is as- 
sumed to be driven by the so-called cosmological pump 
field. However, for our purposes to show that certain 
classes of gravitational waves can induce strong electric 
field perturbations, we use the following solution 



(23) 



h+{t)^e'''^^+^'^exp[bt\ 



with h a constant. Upon substituting into Eq. (I17|) we 
find errors of order 0(l/<^(^^'^) and of the same order 
for Eq. mi). 

In order to show the consistency of the whole model 
we consider Eq. ([TT]) without the first term. This can be 
solved analytically to obtain the exact solution 
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= hotexp 



fcg(l+7)^ 
872 



4t 

i(i+7) 



^l±Jl{t/t,)^'^/('+-'^ 

87 



(24) 



5 



Computing the second derivative we obtain 

1 fc|(i+7) (3,_i)/(i+,) 

+ 27 
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(25) 



Thus the seccond derivative wiU be small with respect 
to the first and zeroth order derivatives, as times passes, 
if the second term, in the first set of brackets, is smaller 
than the first term. This exactly reproduces the condition 
ofEq. m 



B. Large-t Solution 

Here we assume that the last term on the l.h.s of Eq. 
(fT7| is omitted and we have the usual solution of a mode 
that is inside the cosmological Horizon and its amplitude 
decreases with time, 



1 



h+{t) ^ ^ (D^e-^''^' + D^e 



(26) 



When substituted into Eq. (fT7|) we have error terms of 
the order ©(l/t^'+i). Also we find that Eq. ^ is sat- 
isfied to the order 0(1/^4'). 

We can match the two solutions of Eqs. (|24|) and ([26]) 
at t ~ with continuity of the first derivatives and we 
obtain 
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This solution when substituted into Eq. (|13p gives the 
usual oscillating behaviour of the density perturbations 
when the driver (GW) enters the Horizon, after the time 



IV. NUMERICAL RESULTS 



We consider Eq. (flOl) and the first of Eqs. ([9]), where 
we take the magnetic field perturbations equal to zero for 
simplicity. We obtain 
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where the conductivity is in units of (1/fo)- It is evi- 
dent that the conductivity counteracts the electric field 
perturbations that arise from the passage of the GW, 
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FIG. 2: Electric field perturbation as a function of 7 
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whereas the constant h in Eq. ([231) acts in favor of steeper 
in time electric field perturbations. For an initial value 
of {5E^){t[)) = 0.1 we integrate numerically Eq. P5|) for 
CTo = 0-1 J ^0 = 1-0 and 6=1 and study the dependence 
of the electic field perturbation on the parameter 7. This 
is shown in Fig. Also for ctq = 0.1, 6 = 1 and 7 ~ 0.4 
we study the dependence of the electic field perturbation 
on the parameter /iq, effectively the (normalized) am- 
plitude of the metric perturbations. This behaviour is 
shown in Fig. ([3]). 



V. DISCUSSION 

We have considered the propagation of a gravitational 
wave perturbation parallel to the ambient magnetic field 
of an anisotropic magnetized cosmological spacetime, 
non-linearly coupled with its temporal derivatives. It is 
known that in the first order (linear) approximation with 
respect to the amplitude of the perturbation, no elec- 
tric field or density perturbations arise. However in the 
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non-linear level this behaviour changes. Namely while 
certain classes of power-law solutions, for these metric 
perturbations, cannot induce increasing in time electric 
field perturbations, certain other classes of solutions can 
do so. 

Our numerical results and their theoretical counter- 
parts show that while power-law solutions, such as those 
of Eqs. (PTjl and (P^ . are not capable of inducing in- 
creasing in time, electric field perturbations, solutions 
of the form of Eqs. ([23|) and (p4|l can do so. Moreover 
these solutions are consistent with the small-t and large-t 
assumptions, as these are encoded in Eq. p^ . The con- 
sistency of the whole theoretical model is further shown 
by the numerical treatment of Eq. pO|) . The electric 
field perturbations depend on four parameters: First on 
the conductivity ctq, which in an evident and expected 
manner counteracts the generation of electric field per- 
turbations. Second on the constant h, in the class of 
solutions of Eq. (j23p which again in an expected man- 
ner acts in favour of steeper in time electric field and 
density perturbations. The (normalized) amplitude 
of the metric perturbations is the third parameter, as it 
is shown in Fig. ^ and 7 is the fourth parameter, as 
it is shown in Fig. Since the Hubble parameter in 

the z— direction is Hz := (^) = (1 — 7)/(l + 7)i we see 
that higher values of the 7— parameter imply smaller val- 
ues of this parameter, which in turn acts in favour of the 
electric field perturbations. 

Our results hopefully shed some light on the full non- 
linear behaviour of conductive plasmas, in anisotropic 
cosmological models, coupled with an ambient magnetic 
field and driven by (not necessarily weak) metric pertur- 
bations. 
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Also we note that the relativistic Alfven velocity is 

u\ = v\l{l -\-v\) where in conformity with the notation 

used in ^3,] we have 

{F^f (1~7)(37-1) 
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Appendix B 

We present here two of the components of Einstein's 
tensor, just for the purpose of displaying the relative com- 
plexity of the task, computed from Eq. ([T]), after they 
have been carefully expanded up to the second order with 
respect to the GW amplitude. This is a difhcult task and 
has been performed with great care. We obtain 



Go 



and 






2 




hi' 











(33) 



A' 



A 
A 



AW 
AW 



1 



Appendix A 
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We present here for convenience the class of anisotropic 
magnetized cosmologies, with perfect fluid content, 
known as Thome's models Q . The scale factors are given 
by A{t) = a^j^j ^ ^gg^j -^^ ^ jjere 
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The matter energy density and pressure are given by 
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p{t) = 
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Finally the magnetic field points in the 2— direction and 
is given by 
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B,{t)-—, Ho.- (31) 
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Here a dot denotes derivative with respect to the time 
t, while a prime denotes derivative with respect to the 
spatial variable z. Similar expressions are obtained for 
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